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We propose a novel concept for efficient dynamic tuning of optical properties of a high refractive 
index subwavelength nanoparticle with a magnetic Mie-type resonance by means of femtosecond 
laser radiation. This concept is based on ultrafast generation of electron-hole plasma within such 
nanoparticle, drastically changing its transient dielectric permittivity. This allows to manipulate by 
both electric and magnetic nanoparticle responses, resulting in dramatic changes of its extinction 
cross section and scattering diagram. Specifically, we demonstrate the effect of ultrafast switching- 
on a Huygens source in the vicinity of the magnetic dipole resonance. This approach enables to 
design ultrafast and compact optical switchers and modulators based on the ’’ultrafast magnetic 
light” concept. 


Introduction~All-dielectric ’’magnetic light” 
nanophotonics based on nanoparticles of high refractive 
index materials allows manipulation of a magnetic 
component of light at nanoscale without high dissipative 
losses, inherent for metallic nanostructures This 

’’magnetic light” concept has been implemented for 
nanoantennas Q, photonic topological insulators 3, 


broadband perfect reflectors m lll| . waveguides [l3l| . 
cloacking 15, 3; harmonics generation 121, wave-front 
engineering and dispersion control tni. 

Such magnetic optical response originates from the cir¬ 
cular displacement currents excited inside the nanopar¬ 
ticle by incident light. This opens the possibility of in¬ 
terference between magnetic and electric modes inside 
the dielectric nanoparticle at some wavelength. One of 
the most remarkable effects based on this concept is for¬ 
mation of the so-called Huygens source, scattering for¬ 
ward the whole energy [l8| , while for another wavelength 
range, the nanoparticle can scatter incident light in al¬ 
most completely backward direction [13, [ 23 . There¬ 
fore, manipulation by both electric and magnetic reso¬ 
nances paves the way for effective tuning of the dielectric 
nanoparticle scattering in the optical range. The spectral 
positions of the electric and magnetic dipole resonances 
depend on the dielectric particle geometry and ambient 
conditions iiSEEll E2I 1 . Another approach for the 
resonances tuning is to change dielectric permittivity of 
the particle, which was achieved by means of annealing 
of amorphous silicon nanoparticles (^ . 

However, modern technologies require fast, large, and 
reversible modulation of optical response of ultracom¬ 
pact functional structures. For this purpose, different 
types of optical nonlinearities both in metallic 2^ and 
dielectric structures l25| ha ve been utilized such as Kerr- 
e nonlinearities |2fil429l |. free carrie rs g eneration lig¬ 
and variation of their temperature [33| , as well as rel¬ 
atively slow thermal nonlinearity 3^. Since plasmonics 
has high inherent losses and photonic crystals are much 
larger than the wavelength, it is more desirable to use the 
’’magnetic light” concept, dealing with both low-loss and 
subwavelength structures. Moreover, nonlinear manipu- 
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FIG. 1. A schematic illustration of ultrafast scattering ma¬ 
nipulation by a single weak (A) and intense (B) femtosecond 
laser pulse. Intense femtosecond laser pulse switch-on a Huy¬ 
gens source regime, when incident light is scattered in forward 
direction. 


lating by the scattering properties of a nanostructure via 
its magnetic response gives an additional efficient tool for 
ultrafast all-optical switching and light modulation. 

In this work we propose a novel concept of ultrafast 
manipulation of the electric and magnetic responses of a 
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high-index dielectric nanoparticle, employing its ultrafast 
photoexcitation. The existence of the magnetic response 
makes it possible to switch-on the Huygens source (with 
cardioid scattering diagram) in the irradiated nanoparti¬ 
cle on femtosecond time scale (Fig.[T]), enabling to design 
Huygens metasurface with ultrafast laser-induced trans¬ 
parency. In our experimental study, we prove that the 
novel concept ’’ultrafast magnetic light” can be realized 
in non-destructive regime of light-mater interaction. 


Ultrafast silicon photoexcitation-Femtosecond 

(fs) laser pulses are known to provide strong photo- 
induced electronic excitation of free carriers in diverse 
materials (owing to usually low electronic heat capaci¬ 
ties) , which is accompanied by dramatic variation of their 
optical characteristics on the timescale of the laser pump¬ 
ing pulse (< 100 fs). Such almost prompt, fs-laser in¬ 
duced optical tunability appears to be much broader for 
(semi)insulating materials with very minor initial car¬ 
rier concentrations, extending in a sub-ablative regime 
in the visible range, e.g., in terms of dielectric permit¬ 
tivity (e) from large positive (Re(e) ~ -1-10 for the ini¬ 
tial undoped material) to deeply negative (Re(e) ~ —10 
for its strongly-excited, metallized surface lay er) values 
(the so-called insulator-metal transition [35j). In de¬ 
tail, such ultrafast transient modulation of optical di¬ 
electric permittivity in semiconductors and dielectrics is 
related to transient variation of free-carrier (electron-hole 
plasma, EHP) density (peh) through its basic intraband 
and interband contributions 135144,11 1 . Specifically, an op¬ 
tically driven insulator-conductor transition occurs in 
such materials, when the fs-laser fluence-dependent EHP 
frequency passes through the probe frequency, which 
becomes evident as subsequent reflectivity increase for 
stronger ionized materials through their intraband elec¬ 
tronic transitions Simultaneously, ultrafast tran¬ 

sient optical modulation for silicon U, and other 
materials [s^ is additionally enhanced due to a strong 
prompt EHP-driven isotropic renormalization of their di¬ 
rect bandgap, resulting in drastic enhancement of inter¬ 
band transitions and corresponding red spectral shift of 
the optical dielectric permittivity |3,5l | . 


In this work a realistic dependence of the dielectric 
permittivity for photo-excited silicon versus incident fs- 
laser fluence at 800-nm laser wavelength was obtained 
through measurements of single-shot fs-laser pump self¬ 
reflectivity (R) from smooth Si surface at its s- (i?s(45°)) 
and p- (i?p(45°)) polarizations at the 45°-incidence an¬ 
gle and variable effective (absorbed) laser fluences Fgg = 
(1 — i?s,p(45°, E))-F, where F is the incident fluence 
(see Supplementary materials for details). Fitting of 
the experimental dependencies of reflection coefficients 
i?s,p(45°, Eefi) on fluence using a model transient dielec¬ 
tric permittivity for photo-excited silicon. Commonly, it 
can be written as a function of peh as a sum of interband- 
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FIG. 2. (A) Experimental pump self-reflectivity depen¬ 

dences i?s,p(45°, Fefi) (dark squares and circles, respectively) 
with their corresponding fitting red and blue model curves 
I?s,p(45°, Peh). (B) Derived dependences of real (black curve) 
and imaginary (red curve) components of optical dielectric 
permittivity on at the 800-nm pump wavelength. 


and intraband-transition based terms 3^, 42. 431: 


e(a;,peh) = eiB(w*) ( I - 

Wpi(Peh) 


UJ^ + l/{T^{peh)) \ UJTeipeh) 
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Peh 

Phf 
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( 1 ) 


where the prompt peh-dependent bandgap shrinkage ef¬ 
fect on interband transitions is accounted by the spec¬ 
tral dielectric permittivity dependence with the effective 
photon frequency ui* = uj + 0peh/pbgr with the fac¬ 
tor 0, the characteristic renormalization EHP density 
(pbgr), the characteristic band capacity (pbf) of the spe¬ 
cific photo-excited regions of the first Brillouine zone in 
the k-space, the EHP frequency Wpi and electronic damp¬ 
ing time Te at the pump frequency uj (for details of cal¬ 
culations, see Supplementary materials). The resulting 
model peh-dependent pump reflectivities i?s.p(45°, peh) at 
the 800-nm pump wavelength were calculated, using a 
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common Fresnel formula and accounting for all basic ef¬ 
fects (bandgap renormalization^ band filling, EHP screen¬ 
ing of the ion core potential 35, Si), fitting well the 
experimental reflectivity dependences i?s^p(45°,Feff) in 
Fig. [21A. with the characteristic initial dip and the follow¬ 
ing rise. The initial dip implies the predominant intra¬ 
band (Drude) contribution to the optical dielectric per¬ 
mittivity, while the succeeding reflectivity rise indicates 
the presumably interband (red spectral shifting) contri¬ 
bution in the normal dispersion region at the ” red” shoul¬ 
der of the i?i-band of silicon [i^. Such reasonable fitting 
in Fig. [2]A. provides an important relationship between 
magnitudes peh and Fgs in the region, covering the reflec¬ 
tivity dip and rise, which were used to plot the derived 
optical dielectric permittivity versus Fes (Fig. [2l3). Its 
real and imaginary components at the pump frequency 
exhibited versus peh the non-monotonously changing - 
decreasing and then increasing - negative real part and 
monotonously increasing imaginary part (Fig. [2j3), re¬ 
spectively. Such dropping Re(£) and rising Im(£) magni¬ 
tudes are consistent with the increasing intraband contri¬ 
bution at the monotonously rising EHP density, while the 
following rise represents characteristic red-shift changes 
of these quantities across absorption bands in strongly 
photo-excited semiconductors due to their pro mpt peh- 


dependent bandgap renormalization |35l . I41 


Ultrafast photoexcitation of a silicon nanopar- 
ticle-The observed large, ultrafast changes of optical di¬ 
electric permittivity properties in silicon at fluences be¬ 
low its ablation threshold (for spallation under these ex¬ 
perimental conditions, Feff ~ 0.3 J/cm^ 0 ) can sig¬ 
nificantly alter optical response of a silicon nanoparti¬ 
cle, supporting a magnetic Mie-type resonance. Basing 
on the extracted dielectric permittivity values of pho- 
toexcited silicon, we study ultrafast dynamics of scatter¬ 
ing properties (cross section and diagram) of a silicon 
nanoparticle near its magnetic resonance by means of 
full-wave numerical simulations carried out in CST Mi¬ 
crowave Studio. 

We numerically analyzed optical properties of a spher¬ 
ical (the diameter D = 2R = 210 nm) silicon particle 
with its dilectric permittivity, depending on laser fluence 
as shown in Eig. [Tj3. The chosen nanoparticle diameter 
corresponds to excitation of a magnetic dipole Mie-type 
resonance in the vicinity of the femtosecond laser wave¬ 
length A « 800 nm. Its extinction cross-section and scat¬ 
tering diagram are well-known to be rather spect rally 
sensitive in the vicinity of the magnetic resonance [9|, 1201 . 
In particular, at some wavelengths, where magnetic and 
electric dipoles induced in the nanoparticle are almost 
equal and oscillate in phase, the silicon nanoparticle 
works as a H uyg ens source with suppressed backward 
scattering [^, [ij. On the other hand, the scattering 
diagram can be tuned to the regime of suppressed for¬ 
ward scattering (’’reverse” Huygens source), when mag¬ 
netic and electric dipoles oscillate with phase difference 
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FIG. 3. Numerically calculated extinction spectrum of a sili¬ 
con sphere of a radius D = 210 nm (A). Insets: numerically 
calculated electric and magnetic fields distributions near the 
silicon sphere. Scattering diagrams of the silicon sphere at A 
= 740 nm (B), 800 nm (C) and 860 nm (D). Black arrows in 
(B-D) indicate direction of light incidence. 


of 7r/2 @,[ 23 . 

The results of numerical simulations of extinction and 
scattering properties of the silicon nanoparticle in the 
vicinity of the magnetic dipole resonance are shown in 
Fig. El The calculated distributions of electric and mag¬ 
netic fields prove emerging a magnetic dipole moment 
in this range (Eig. EK)- Since both extinction spec¬ 
tra (Eig. EK) and scattering diagram (Eig. EP-D) are 
strongly wavelength-dependent, it is possible to tune 
these two parameters by varying either the nanoparticle 
size or its dielectric permittivity. 

In order to considerably tune silicon nanoparticle op¬ 
tical properties near the magnetic dipole resonance, it is 
necessarily to have time of the mode formation shorter 
than the laser pulse duration. For D = 210 nm sili¬ 
con nanoparticle, the magnetic resonance mode has ap¬ 
proximately 50-nm the full-width at half-maximum, and, 
therefore, the Q-factor is about 16, corresponding to the 
mode formation time in such open resonator of about 7 
fs. Since the typical laser pulse duration of the most com¬ 
mercial femtosecond systems is more than 30 fs and op¬ 
tical changes are stabilized at the middle of the fs-pulse, 
the character of interaction between the nanoparticle and 
the most part of the femtosecond pulse is significantly 
changed in comparison with the case of nanoparticle ir- 
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FIG. 4. Numerically calculated normalized extinction spectra (A) and scattering diagram (B) of the 210-nm large silicon 
particle irradiated at A = 800 nm and FeS = 0 (black curve), 35 mj/cm^ (red curve), and 100 mj/cm^ (blue curve). (C) 
Numerically calculated dependences of reflectance (black curve) and transmittance (red curve) through array of such 210-nm 
silicon spheres with the period of 600 nm on effective fluence F’eff. 


radiation by low intensity. Therefore, this concept can 
be applied for effective ultrafast nanophotonic devices 
(switchers, modulators, etc.). 

To simulate numerically changes in optical properties 
of photoexcited silicon sphere with the diameter of D 
= 210 nm at the wavelength A = 800 nm, we use the 
derived above dependencies of Re(e) and Im(e) on ab¬ 
sorbed laser fluence Tea (Fig. 03). The considered range 
of absorbed fluences Fgg < 100 mJ/cm^ corresponds to 
non-ablative regime of the laser-particle interaction, but 
laser fluences are still high enough to generate rather 
dense EHP (peh~lxl0^^ cm“^) for efficient switching of 
the optical properties of the nanoparticle. As was men¬ 
tioned before, the scattering diagram of the 210-nm large 
nanoparticle is almost symmetric at A = 800 nm and 
Fes ~ 0 (Fig. EP), while the extinction cross section 
((Text) normalized on ttR^ has rather high value of about 
9 (see Figs.EP and HP). The latter parameter is changed 
almost by three times with fluence increasing up to Tgff 
= 100 mJ/cm^ at the fixed wavelength 800 nm, owing 
to the strong shift of the peak position of the extinction 
spectrum (Fig.EJA). Its scattering diagram appears to be 
also very sensitive to the corresponding changes of the 
dielectric permittivity. For D = 210 nm at A ~ 810 nm, 
the transition from the typical dipole scattering diagram 
to the Huygens source forward scattering is observed in 
the fluence range Fgs = 0 — 100 mJ/cm^ (Fig.|3|3). 

Since the scattering cross section is reduced and back¬ 
ward scattering is suppressed, almost all energy of the 
incident light is transmitted through the nanoparticle at 
negligible reflection. It means that a metasurface based 
on such elements should demonstrate ultrafast switch¬ 
ing from a reflecting to a non-reflecting state. In order 
to propose the design for the ultrafast switching meta¬ 
surface, transmission and reflection are numerically cal¬ 
culated for a infinite periodic array of silicon spheres 
with the diameter H = 210 nm and period of 600 nm. 
The resulting dependences R(Feff) and T{Fes) for the 


800-nm wavelength demonstrate significant changes even 
at Eeff = 5 mJ/cm^ where transmission is increased by 
more than 100% (Fig. |4p). Taking into account 10-fold 
average intensity enhancement within each nanoparticle 
(Fig. EP), the estimated incident fluence, corresponding 
for 100% increasing of transmission, is about F Ri 0.5 
mJ/cm^. For excitation using a diffraction-limited spot 
of around 1 area, the above fluence gives a switching 
energy of 5 pJ, being comparable with microphotonic ring 
resonators (25 pJ) l30| and different plasmonic nanos¬ 
tructures (7 - 20 pJ) Moreover, despite the in¬ 

crease of Im(e) with fluence (Fig. EB); the transmission 
of the silicon nanoparticle-based metasurface is increased 
by more than 1000% at Fes = 100 mJ/cm^ (Fig.EP) due 
to switching into the Huygens source scattering regime. 
Therefore, the proposed concept of ’’ultrafast magnetic 
light”, dealing with both extinction cross section and 
scattering diagram control, provides low-intensity non¬ 
linear light manipulation by means of ultimately simple 
nanoobject - a silicon nanosphere. 

Conclusion-In summary, ultrafast generation of 
dense electron-hole plasma in a dielectric nanoparticle, 
supporting a magnetic dipole resonance in the optical 
range, paves the way for the novel concept of ’’ultrafast 
magnetic light” based not only on tuning of extinction 
cross section of the nanoparticle, but also on tuning of 
its scattering diagram. In the frame of this concept, the 
possibility of ultrafast switching on a Huygens source in a 
silicon nanoparticle photoexcited by a femtosecond laser 
pulse is shown, enabling high-efficient ultrafast light ma¬ 
nipulation. Tuning its dielectric permitivity via fs-laser 
photoexcitation will open a novel class of ultrafast de¬ 
vices, based on the diversity of ’’magnetic light” effects. 
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Supplementary materials for ’’Ultrafast Mag¬ 
netic Light” 

Experiment on ultrafast silicon photoexcitation 

In the experiments we used a commercial 0.45-mm 

thick, atomically flat undoped silicon Si(lOO) wafer with 

a few nanometer-thick native oxide layer, arranged on a 
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three-dimensional motorized translation micro-stage un¬ 
der PC control and raster-moved from laser shot to shot 
to expose its fresh surface spots. Its single-shot (one 
laser pulse per surface spot) laser irradiation at a repe¬ 
tition rate of 10 Hz was provided by single IR (800 nm) 
Ti: sapphire laser pulses in TEMqo mode with a FWHM 
(full-width at half-maximum) duration Tp Ri 100 fs and 
pulse energies up to 1.5 mJ, using a triggered electro¬ 
mechanical shutter. The pump pulse energy was re¬ 
duced in these experiments, using a combination of a 
half-wave plate and a Glan-prism polarizer, to the min¬ 
imum level slightly above i?max ~ 0.2 mJ to avoid the 
laser beam degradation due to its self-focusing in air and 
air plasma scattering/refraction (the critical power at the 
wavelength « 3 GW [i^, i.e., « 0.3 mJ for the 100-fs 
laser pulses). 

The fs-laser pump pulses were focused in s- and p- 
polarizations at the angle of 45° by a silica lens (focal 
distance f=500 mm) into a small focal spot (the Gaus¬ 
sian main ellipse 1/e-radii ~ 0i/e,x ~ 0-10 mm and 
0i/e,y ~ 0.05 mm) on the surface of the Si sample. 
Energies of of incident pump laser pulses Ei < i?„iax 
were monitored by means of a pyroelectric energy meter 
(OPHIR), respectively, while energies of correspond¬ 
ing reflected pump laser pulses were simultaneously mea¬ 
sured 3-5 times for each pulse energy over fresh surface 
spots by another pyroelectric energy meter (OPHIR), 
providing the pulse energy dependence of the pump self¬ 
reflectance i?s,p(45°, Ei). 

In the experimental arrangement of incident and re¬ 
flected energy measurements (Ei and Es,p, respectively), 
the measured self-reflectivity values i?s,p(45°, Ej) repre¬ 


sent their true i?s,p(45°, li), temporally averaged over the 
100-fs wide pump profile /(t) and over the laser spot laser 
spot with the spatial distribution /(x,y) characterized by 
the peak magnitude li 


i?,.p(45°,Ei) = 


_ III Rs,’p{'i^^°,h)h{x,y,t)dxdydt 


Ei 


/// Ii{x,y,t)dxdydt 


( 2 ) 

where x,y are the spatial surface coordinates. The tempo¬ 
ral averaging hides the unknown dynamics of EHP in the 
photo-excited silicon, while the latter spatial averaging 
can readily be lifted up via a common iterative decon¬ 
volution procedure 35, 4^, employing Es_p(45°,E^) 

values measured at lower energies Ei (peak fluences Ej) 

0. 


E,,p(45°,E0 = 

Es.p(45°,E0 - E i?,,p(45°,Efe)Efe(E„Efc,Efe_i) 

__ ('3') 

E,(Ei,Ei,Ei_i) ’ 

where local reflectivity values across the laser spot are 
assumed to be constant in each range Ek-i ^ Efc and are 
weighted with their statistic weights Pk{Ei, Fk, Ek-i), 
which are equal at each pulse energy Ei to the partial 
energies in the ring-like squares ASk{Ei, Fk, Fk-i) with 
the boundary local fluences Fk-i and Fk 

Tl T7 T7 ^ AEfc(Ej,Efe,Efc_l) FkASk(Ei,Fk,Fk-l) 

Efc(E„Efc,Efe_i) =- — ---- 

jZ FkASk{E,,Fk,Fk-i) 

k—1 1 /A 



FIG. 5. Experimental setup: BS-beam splitter, AC- 
autocorrelator, QPW-|-GPP-energy attenuating combination 
of a quarter-wave plate and a Glan-prism polarizer, EM 1,2- 
thermocouple energy meters, M-mirror, L-focusing silica 
lenses, DO-digital oscilloscope, 3D-MS-three-dimensional 
motorized micro-stage, PC-laptop for data acquisition and 
hardware control. 


The spatial deconvolution procedure starts from the 
lowest experimental Ei (Ei) magnitude, where the depen¬ 
dences i?s.p(45°, Ej) exhibit their almost unperturbed re¬ 
flectivity values [4J, and proceeds eventually to higher 
fluences. In both the cases, the incident fluence cali¬ 
bration, providing the peak surface fluence F up to 1.5 
J/cm^ as a function of the incident energy and the above- 
mentioned Gaussian beam parameters, was performed by 
means of an optical microscope to measure main radii of 
the resulting ablative surface craters. 

Details on modeling of ultrafast silicon photoex¬ 
citation 

During modeling of ultrafast transient optics of pho- 
toexcited silicon with the help of Eq. 1 in the main 
text, important instantaneous electronic bandgap renor¬ 
malization and screening effects were taken into account 
for the first time. Also, we used the following param¬ 
eters: the characteristic renormalization EHP density 
Pbgr ~ 1 X 10^^ cm“^ 0, being typically about 5% of 
the total valence electron density (r; 2 x 10^^ cm“^ in Si) 
to provide the ultimate 50% electronic direct bandgap 
renormalization T^, i.e., TiQ k, 1.7 eV of the effective 
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minimal gap ~ 3.4 eV in silicon l48|, while pbf is 
the characteristic band capacity of the specific photo- 
excited regions of the first Brillouine zone in the k-space 
(e.g., pu{L) Ri 4 X 10^^ cm“^ for L-valleys and Pbf(^) ~ 
4.5 X 10^^ cm“^ for X-valleys in Si), affecting interband 


transitions via the band-filling effect 
bulk EHP frequency Wpi is defined as 


Mka. The 


Wpi(peh) 


_Pehe^_ 

EO^hf (Peh)nT.opt (peh) 


(5) 


where the effective optical (e-h pair) inass m* ^ r; 0.14me 
in L-valleys or 0.19 in X-valleys [s^, IH, |3 SS SI] is 
a peh-dependent quantity, varying versus transient band 
filling due to the band dispersion and versus bandgap 
renormalization (H^ . The high-frequency electronic di¬ 
electric constant Ehf was modeled in the form ehf(peh) = 
1 + Ehf (0) X exp(—peh/pscr), where the screening density 
Pscr ~ 1 X 10^^ cm“^ was chosen to provide ehf —>■ 1 in 
dense EHP. The electronic damping time Te in dense EHP 
at the probe frequency Wpr was taken, similarly to met¬ 


als, in the random phase approximation as proportional 
to the inverse bulk EHP frequency 


Te 


f l2SEl 

yTT^ySWpl 


1 -I- exp 


hijJ 

ksTe 


(nksTe)'^ -I- {hujy 


( 6 ) 


where Ep Ri 1 - 2 eV is the effective Eermi-level quasi¬ 
energy for electrons and holes at peh < 1 x 10^^ cm“^, 
h and ks are the reduced Planck and Boltzmann con¬ 
stants, respectively, and Te is the unified EHP tempera¬ 
ture, being a weak function of peh [l^- Here, the latter 
relationship was evaluated for Tllo > fceTe in the form 
Te(peh) ~ 3 X 10^/(a;pi(peh)w^), accounting multiple car¬ 
rier scattering paths for the three top valence sub-bands, 
and multiple X-valleys in the lowest conduction band of 
silicon. 

The resulting peh-dependent model oblique-incidence 
pump reflectivities i?s_p(45°,Peh) calculated, using com¬ 
mon Fresnel formulae, fit well the extracted experimental 
reflectivity dependences i?s,p(45°, Teff) in Fig. 2A with 
the characteristic initial dip and the following rise. 












